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Prime factorization is an outstanding problem in arithmetic, with important consequences in a 
variety of fields, most notably cryptography. Here we employ the intriguing analogy between prime 
factorization and optical interferometry in order to obtain, for the first time, analytic expressions 
for closely related functions, including the number of distinct prime factors. 
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For millennia, people have been fascinated with the 
world of numbers and in particular with prime numbers. 
As early as 300 BC Euclid proved the infinitude of the 
primes and laid the foundation for the proof of the fun- 
damental theorem of arithmetic. The last centuries have 
seen a surge in scientific activity focused on the primes 
and related mathematical themes, with leading mathe- 
matician such as Euler and Riemann greatly contribut- 
ing to our current knowledge regarding the complexity 
of the distribution of primes. However, prime factoriza- 
tion of integers, especially that of large numbers, is still 
a formidable task. As a consequence, different applica- 
tions in fields such as cryptography have been developed, 
which take advantage of the great difficulty of factoring 
large numbers (e.g., the RSA cryptosystem [T]). 

During the last decades, physicists have contributed 
to the investigation of integer factorization as well as re- 
lated mathematical topics, such as the Riemann hypoth- 
esis (see [2] and the references therein), through investi- 
gating the relationship with various physical systems. In 
particular, the intriguing relationship between quadratic 
Gauss sums and integer factorization has led to recent 
experimental realization of integer factorization through 
NMR and optical interference (see e.g., [SHZ]). 

In this Letter we explore the relationship between 
prime factorization and optical intereferometry with the 
aim of obtaining novel analytic expressions for number- 
theoretic functions directly related to prime factoriza- 
tion. Our study, which is based on the relationship be- 
tween the multiple-slit interference experiment and the 
occurrence of the primes, exhibits the potential of physi- 
cal analogy not only in realizing a mathematical enigma 
experimentally, but also contributing to the theoretical 
endeavor to unravel the enigma. 



The fundamental theorem of arithmetic states that any 
positive integer n can be decomposed into prime numbers 
in a unique manner, that is: 



n 


2 


3 


5 


7 


11 


13 


17 


19 


23 




1 
































2 


1 


























1 


3 





1 























1 


4 


2 


























1 


5 








1 




















1 


6 


1 


1 























2 


7 











1 

















1 


8 


3 


























1 


9 





2 























1 


10 


1 





1 




















2 


11 














1 














1 


12 


2 


1 























2 


13 

















1 











1 


14 


1 








1 

















2 


15 





1 


1 




















2 


16 


4 


























1 


17 




















1 








1 


18 


1 


2 























2 


19 























1 





1 


20 


2 





1 




















2 


21 





1 





1 

















2 


22 


1 











1 














2 


23 


























1 


1 


24 


3 


1 























2 


25 








2 




















1 



TABLE L Prime factorization of the first twenty-five positive 
integers. The rightmost column presents the corresponding 
number of distinct prime factors. 



uj(n) 

Pi , (1) 

i=l 

where cij(n) is the number of distinct primes function, pi is 
the the i*^ prime and ap. (n) is the corresponding power. 
Table |I] depicts the integer decomposition into primes of 
the first twenty- five integers. From the table one can 
deduce three features of prime factorization. First, one 
notes that each prime appears periodically as n increases, 
with a period equal to the corresponding prime. This is 
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a consequence of the fact that for a given prime p all 
integers divisible by p can be written as Np = mp, where 
m is an integer. We define the following function, which 
conveys the periodicity of the occurrence of a prime p: 



■4>p{n) = 



1 ap{n) > 
otherwise 



(2) 



The second feature concerns with the synchronization in 
the occurrence of the different primes in Tablejlj or equiv- 
alently, the synchronization of the ipp{nys. A close exam- 
ination shows that the occurrence is in phase with respect 
to extrapolating the periodic set {ipp{n)} backward to- 
ward zero, where "0^(0) — 1 for all primes. The third 
feature we note is that a particular prime number p* ap- 
pears in the table when ipp{0) — for all p < p* . Thus, 
knowledge of an initial set of m primes 2, 3, ■■■,Prn allows 
for the determination of at least one consecutive prime. 
In fact, it is easy to be convinced that this set allows 
for the determination of all consecutive primes obeying: 

P < pi- 
Examining Eq. [2] we recall a similar functional de- 
pendency which appears in optical intcrfcrometry. It is 
well known that the normalized intensity of p coherent, 
monochromatic sources in the form of infinitely narrow 
slits positioned on a line with spacing d is given by [5], 



IpiO) = 



1 sin'^ {p'Kdsm{6)/X) 
p^ sin"^ {ndsm{6)/ X) 



(3) 



where, A and 6 are the wavelength and angle with the 
normal at the midpoint of the slit arrangement, respec- 
tively. Figure [l] depicts the intensity for the particular 
case p = 5, where we have used: x = dsin{9) / X. The 
function Ip{x) can be made to coincide for all positive 
integer values with ^pp{n) through the coordinate trans- 
formation x' = xp. This is analogous to scaling the dis- 
tance between sources by p in Eq. [3] a scaling which will 
prove valuable further on when we superimpose contri- 
butions of different (non-coherent) sets of sources. Thus, 
Ip{x') has its maximal value of unity for any integer di- 
visible by p, and vanishes for all other integers. In what 
follows we will omit the superscript '. 

The fact that 'ipp{n) = Ip{n) enables us to make the 
following statement with respect to the third feature 
pointed out earlier. Let 2,3, ...,Pm be a set of m initial 
primes and define: 



1=1 



Ip, {x) 



m 



sm^ (ttx) 



^ pf sin^ [i^x/pi) 



(4) 



Then finding the consecutive set of primes which obey 
p < p^ is equivalent to finding the zeros of uJmix) in the 
interval Pm < x < p^. This is a direct consequence of the 
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FIG. 1: Normalized intensity corersponding to a set of five 
monochromatic, coherent narrow slits. 



fact Ip{x) obtains its zeros when the numerator vanishes 
but the denominator does not. Thus, Ip{x) = for any 
integer not divisible by p, and therefore Wm(a;) obtains 
its zeros for any integer which is not divisible by the first 
m primes. Figure [5] depicts the case m = 3 (pa = 5), for 
which the zeros of u)3{x) in the interval 5 < x < 25 are 
located at 7, 11, 13, 17, 19 and 23. 




FIG. 2: Functional dependence of 0)3 (x) (pa = 5), which rep- 
resents the intensity corresponding to the set of sources de- 
picted in Fig. [3] The zeros in the interval 5 < a; < 25 are 
located at the consecutive primes: 7, 11, 13 17, 19 and 23. 



The analogy with interference of narrow slit sources is 
conceptually rather simple, though experimental imple- 
mentation might prove difficult. Equation [4] represents a 
superposition of contributions due to m sets of infinitely 
narrow slit sources (in an experimental realization one 
would need to take into consideration the effect of diffrac- 
tion due to the finite width of the slits) . Each such set has 
Pi coherent sources with spacing d/pi. While all sources 
(in all sets) have the same wavelength A, different sets 
are non-coherent. Figure |3] depicts the arrangement of 
narrow slit sources corresponding to the to = 3 scenario. 
Note that due to the fact that the superposition is of 
sets consisting of a prime number of sources, the sources 
never overlap, except for the middle point (marked white 
in Fig. [3]), which is shared by all odd sets. We also note 
that if the different sets (i.e., pertaining to different prime 
numbers) were coherent, which would be equivalent to il- 
luminating the slit arrangement with a single monochro- 
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matic, coherent beam, the functional dependency of the 
intensity would be qualitatively different, in large part 
due to the coupling terms representing interference be- 
tween the different sets. 
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FIG. 3: Arrangement of narrow slit sources (represented by 
circles) on a line segment of length d. The sources are divided 
into three sets of 2 (red), 3 (green) and 5 (blue) sources with 
spacings d/2, d/3 and d/5, respectively. The only overlapping 
point is the midpoint, shared by all odd sets (white circle). 

The partial sum u}m{x) has three qualitatively different 
regions. For 2 < x < Pm, 0Jm{n) faithfully represents 
the number of distinct primes in the factorization of n. 
In particular, ujm{n) — 1 for n = pi. In the second 
region, pm < x < the zeros of u!rn{x) coincide with 
the consecutive set of primes {pm < P < Pm)- the 
third region p^ < x, ojm{n) in general neither represents 
the number of distinct prime factors, nor has its zeros at 
locations of prime numbers. 

There is an intriguing connection between w„i(n) and 
the sieve of Erastosthenes, which provides a rigorous pro- 
cedure for finding all primes smaller or equal to a given in- 
teger N . Starting from the first prime pi — 2, one marks 
all multiples of the prime up to N . The first integer left 
unmarked at every step is by definition a prime number 
and is used in the next step in order to locate additional 
composite numbers. The procedure terminates when one 
has located all primes smaller or equal to ^/N . In a simi- 
lar manner, as pointed out above, the partial sum Wm(n) 
distinguishes between composite integers {uj^in) > 0) 
and primes (w„i("-) = 0) in the interval p„i < n < 
based on a known initial set of m primes. Thus, by lo- 
cating the zeros of a;,„(7i) in this interval one obtains the 
consecutive primes Pm+i, ■■■,Pm+k (where Pm+k is the 
closest prime smaller than The next step in the 

procedure is to locate the zeros of uJm+k {n) in the inter- 
val {pm+k,Pm+k)- sieve of Erastosthenes, the 
procedure ends when one has located the largest prime, 
which is smaller or equal to Vn. 

An important consequence of Eq. [4] is: 



io{n) = lim ujm{n) 
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1 sm^ (im) 



^ sm^ {-Kn/pi) 



(5) 



It is important to note that, apart from being important 
in its own right (e.g., [HI HH]), w(n) is related to other 
important functions encountered in number theory (and 



elsewhere), such as the prime-counting function tt^x) and 
Riemann's zeta function Q{s) (s > 1). In particular, one 
has the following identity [TT] . 
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The analogy between prime factorization and optical 
interferometry can be taken even further. An important 
function related to prime factorization is the total num- 
ber of prime factors (with multiplicity) defined as, 



Sl(ri) = 2, cupi ('^)- 



(6) 



In order to obtain an expression for 51 (n), we are led to 
examine the dependence of ctp over the positive integers 
which are divisible by p (for all other integers ctp = 0). 
We have already seen that Ip{n) conveys the periodic- 
ity of occurrence of the prime p in the factorization of 
integers. 
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FIG. 4: Relationship between Qp(n) and the set {Ip^{Ti)} 
(m=l,2,3,...) for the particular case p = 3. 

Figure |4] depicts the first fifteen non- vanishing entries 
for the particular case Q;3(n), together with the corre- 
sponding values of /s, /g, l27 and /gi. A close examina- 
tion leads to the conclusion that Up is a superposition of 
contributions from: /„, /, 



00 00 
ap(n)=^/p,(n) = ^— - 



siTn? (nn) 



p^i sirfl (nn/p^) 
Substituting Eq. [7] into Eq. |6]one then gets: 

i=i j=i i=i j=i Pi sm^ [T^n/pl 1 



(7) 



(8) 



In summary, we have explored the analogy between 
prime factorization of integers and optical interferome- 
try with the aim of obtaining analytic expressions for 
number-theoretic functions closely related to prime fac- 
torization. It was shown that the normalized intensity 
corresponding to a set of p coherent sources possesses 
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the functional dependency pertaining to the occurrence 
of the corresponding prime in the factorization of the 
integers. We then examined the partial sum ujmin)^ per- 
taining to the intensity field of m sets, each consisting 
of a prime number of narrow slits, and elaborated on its 
connection with the sieve of Erastosthenes. An analytic 
expression for the number of distinct prime factors w(n) 
was then derived through taking the limit of infinite sets 
of slits. Finally, we also obtained analytic expressions for 
the functions ap{n) and fl{n), which together with uj{n) 
completely define prime factorization. 
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